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I. INTRODUCTION 



When an open quantum system undergoes a dynamical 
evolution generated by a Hamiltonian, how far does its 
state evolve from itself as a function of the magnitude of 
the Hamiltonian? This is a fundamental question that is 
central to quantum information science [l| and quantum 
control 0, [3]. In order to make it more precise, suppose 
that the unitary propagator of the evolution, U(t), that 
is related to the Hamiltonian H(t) via the Schrodinger 
equation U = —j^HU, is written in terms of an effec- 
tive Hamiltonian £l(t) as U(t) = exp[— itCl(t)/H]. The 
effective Hamiltonian can in turn by calculated from the 
Hamiltonian H(t), using a Dyson or Magnus expansion. 
If the system S is "open", it is coupled to a bath B, 
and its time-evolved state is given via thepartial trace 
operation by p s (t) = Tr B U(t)p(0)U(t)^ 0, H|. In the 
context of quantum information, the question we have 
stated pertains to the problem of "quantum memory", 
i.e., what is — As(0)|| as a function of ||f2(£)|| or 

\\H (i)||? When the issue is quantum computation or gen- 
eral quantum control, one is interested in comparing two 
time-evolved states: the "ideal" state Pg(t) that is error- 
free and is described by an "ideal" effective Hamiltonian 
(e.g., for a quantum algorithm), and the "actual" 
state ps(t), that underwent the full noisy dynamics de- 
scribed by the total effective Hamiltonian f2(i). Then the 
question becomes: what is ||/0s(£) — Ps(OII as a function of 



\\£l{t)-n°(t)\\ or \\H(t)-H°(t)\\? The memory question 
can of course be seen as a special case of the computation 
question. 

Here we prove bounds that answer these fundamen- 
tal questions. Our bounds have immediate applications 
to problems in decoherence control @ and fault-tolerant 
quantum computation Q, as they quantify the sense in 



which a distance between (effective) Hamiltonians de- 
scribing the evolution should be made small, in order to 
guarantee a small distance between a desired and actual 
state. 

To begin, we first recall the definition and key proper- 
ties of so-called unitarily invariant norms, as we use such 
norms extensively (Section [TTJ) . We mention the trace 
norm and operator norm, and introduce a new norm that 
mixes them. We then briefly review the accepted distance 
measures between states, so as to quantify the meaning of 
an expression such as || ps(t) — P%(t) || (Section [Hi]). Next 
we discuss how to compute the effective Hamiltonian Q,(t) 
using the Magnus expansion or Thompson's theorem, and 
introduce a generalized effective superoperator generator 
( Section IIV[) . Since in many applications one is inter- 
ested not in the distance between states generated by 
unitary, closed-system evolution, but instead in the dis- 
tance between states of systems undergoing non-unitary, 
open-system dynamics, we prove an upper bound on the 
distance between such system-only states, in terms of the 
distance between the full "system plus bath" states (Sec- 
tion [V}. We then come to our main result: an upper 
bound on the distance between system states in terms 
of the distance between (effective) Hamiltonians describ- 
ing the system+bath dynamics (Section |VT|. We present 
a discussion of our result in terms of an example bor- 
rowed from decoherence control using dynamical decou- 
pling fSection IVII[) . We conclude in Section FVIIII with 
some open questions. 



II. UNITARILY INVARIANT NORMS 

Unitarily invariant norms are norms satisfying, for all 
unitary U,V @: 
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||EMV|| ui = ||A|| ui . (1) 

We list some important examples, 
(i) The trace norm: 

P|| lS TrL4|=5>(A), (2) 



where \A\ = \J A^A, and Si(A) are the singular values 
(eigenvalues of \A\). 

(ii) The operator norm: Let V an inner prod- 
uct space equipped with the Euclidean norm ||x|| = 

VJ2i \ x i\ 2 ( e i, e i): where x = J2i x i e i e v and v = 
Sp{ei}. Let A : V V. The operator norm is 



HAJloo,! = sup Vf^= su P W A Ph, (11) 
P ec(H) \\P\\i \\p\\x=i 

where A : V h V, and V is a normed space equipped 
with the trace norm. Note that if Ap is another normal- 
ized quantum state then HAH^i = 1. Also, by defini- 
tion ||Ap||i < 1 1 A|| oo ; i . Moreover, it follows immediately 
from the unitary invariance of the trace norm that the 
O-T norm is unitarily invariant. Indeed, let T = V ■ 
be a unitary superoperator (i.e., V is unitary); then 

HrxArJlU,! = supJ^Ap^Mlli = supJIApllx = 
||A||oo,i- Therefore the O-T norm is also submultiplica- 
tive. However, note that unlike the case of the standard 
operator norm, there is no simple expression for HAHo^i 
in terms of the singular values of A. 



|| Alloc = sup = max Si (A). (3) 

xev IfII * 

Therefore ||Ax|| < 1 1 A 1 1 1 1 a; 1 1 . In our applications V = 
C(TL) is the space of all linear operators on the Hilbert 
space 7i, A is a superoperator, and x = p is a normalized 
quantum state: ||p||i = Trp = 1. 

(iii) The Frobenius, or Hilbert-Schmidt norm: 



\\A\\ 2 = ^m= ^ Si {Af. (4) 

All unitarily invariant norms satisfy the important prop- 
erty of submultiplicativity Q: 

\\AB\U < \\A\U\B\U. (5) 

The norms of interest to us are also multiplicative over 
tensor products and obey an ordering 0: 

\\A®B\\i = \\AU\B\U i = l,2,oo, 
Plloo < ||A|| 2 < HAIU, 
\\AB\\ ui < HAHoollflllui.llBlUIIAIU. (6) 

There is an interesting duality between the trace and 
operator norm [9J: 

PHi = max{\Tr(B*A)\ : < 1} (7) 

HAIloo = max{|Tr(BtA)| : ||B|| a < 1}, (8) 
\Tr(BA)\ < IIAIIoollBtll^llBtll^UAIl! (9) 

In the last three inequalities A and B can map between 
spaces of different dimensions. If they map between 
spaces of the same dimension then 

P||i = max |Tr(BU)|. (10) 

B—I 

We now define another norm, which we call the 
"operator-trace" norm (O-T norm): 



III. DISTANCE AND FIDELITY MEASURES 

Various measures are known that quantify the notion 
of distance and fidelity between states. For example, the 
distance measure between quantum states p\ and p2 is 
the trace distance: 

D{px, P2 )= gHpi-palli. (12) 

The trace distance is the maximum proba- 
bility of distinguishing pi from p 2 - Namely, 
D{pi,p 2 ) = max <p</((P)i - (P)2), where 
(P)i = TrpiP and P is a projector, or more gener- 
ally an element of a POVM (positive operator- valued 
measure) [l|. The fidelity between quantum states p\ 
and p2 is 



F(pi,p 2 ) = \\y/p7-s/p2~\\l = \l yJ~Plp2y/~Pl, (13) 

which reduces for pure states p\ = \ip) {ip\ and 
P2 = |0) (0| to F(\i/>) , |0)) = | (VI <t>)\- The fidelity and 
distance bound each other from above and below [10| : 

1 - D{p 1 ,p 2 ) < F(pi, to) < Vl - D{ Pl , p 2 )\ (14) 

so that knowing one bounds the other. Many other mea- 
sures exist and are useful in a variety of circumstances 

a. 

IV. GENERATORS OF THE DYNAMICS 

A. Effective superopator generators 

We shall describe the evolution in terms of an effective 
superoperator generator L{t), such that 

p(0) ^ p(t) ee e tL Vp(0). (15) 
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The advantage of this general formulation is that it in- 
corporates non- unitary evolution as well. Nevertheless, 
here we focus primarily on the case of unitary evolution 
p(t) = U{t)p{Q)U{t)\ with U = -%HU, for which we 
have 



L(t) 



:[«(*).• 



(16) 



This follows immediately from the identity e lA pe lA = 
e ~ l l A '-]p = J2^=o ^~n\ [nA, p]> satisfied for any Hermitian 
operator A, where [ n A, p] denotes a nested commutator, 
i.e., [ n A,p] = [A, [ n -iA, p]], with [ A,p] = p. 



B. Magnus expansion 

In perturbation theory the effective Hamiltonian can 
be evaluated most conveniently by using the Magnus ex- 
pansion, which provides a unitary perturbation theory, in 
contrast to the Dyson series [ll[ . The Magnus expansion 
expresses f2(t) as an infinite series: Q(t) = J2^Li ^n(t), 
where Oi(t) — j f H{t\)dt\, and the nth order term in- 
volves an integral over an nth level nested commutator 
of H(t) with itself at different times. A sufficient (but 
not necessary) condition for absolute convergence of the 
Magnus series for Q(t) in the interval [0, t) is [1 21 ] : 



\H(s) 



, ds < 7T. 



(17) 



C. Relating the effective Hamiltonian to the "real" 
Hamiltonian 



There is also a way to relate the effective Hamiltonian 
to the real Hamiltonian in a non-perturbative manner. 
To this end we make use of a recently proven theorem 
due to Thompson [H,[l|]. 

Let us consider a general quantum evolution generated 
by a time-dependent Hamiltonian H , where V is consid- 
ered a perturbation to Hq (in spite of this we will not 
treat V perturbatively): 



H(t) = H (t) + V(t). 
The propagators satisfy: 
dU(t,0) 



dt 



dU o (t,0) 
dt 



-iH(t)U{t,0), 



-iH (t)U (t,0). 



(18) 



(19) 



(20) 



We define the interaction picture propagator with respect 
to Hq, as usual, via: 



U(t,0) = U (t,0^U(t,0). 



(21) 



It satisfies the Schrddinger equation 
dU(t,0) 



dt 



iH(t)U(t, 0), 



with the interaction picture Hamiltonian 
H(s) = U (t,0^V(t)U (t,0). 



(22) 



(23) 



See Appendix [A] for a proof. 

The interaction picture propagator U (t, 0) can be for- 
mally expressed as 



T 



exp 



H(s)ds 



exp[— itCl(t)], 



(24) 
(25) 



where the second equality serves to define the effective 
interaction picture Hamiltonian Q(t). 



Lemma 1 There exist unitaries {W(s)} such that 



f2(i) = - / W(s)H(s)W(s)Us. 



(26) 



This is remarkable since it shows that the time- 
ordering problem can be converted into the problem of 
finding the (continuously parametrized) set of unitaries 
{W(s)}. 

Proof. The formal solution (|2"4")l can be written as an 
infinite ordered product: 



U(t, 0) = lim TT exp 



3=0 



N \N 



(27) 



Thompson's theorem [l~3|, [l4[ states that for any pair of 
operators A and B, there exist unitaries V and W, such 
that e A e B = e VAV^+WBw\ It fou ows immediately that 
if {Aj}jL are Hermitian operators then it is always pos- 
sible to find unitary operators {W^ jjLi such that 

N N 

|Q expl-iAj] = expl-i^WjAjW}]. (28) 

The proof is non-constructive, i.e., the unitaries {Wj}^ =l 
are not known. Applying this to Eq. (|27[) yields 



U(t, 0) = ^lim exp H E W,B [§)wj] 



= exp[-i / W(s)H(s)W(s) f ds}, (29) 
Jo 

which is the claimed result with the effective Hamiltonian 
Cl(t) defined as in Eq. ((26]). ■ 

An immediate corollary of Lemma [T] is the following: 
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Corollary 1 The effective Hamiltonian Q(t) defined in 
Eq. \2b\) satisfies, for any unitarily invariant norm: 

< ~Jd8\\V(a)U = (\\VU (30) 



< sup \\V(s)\\ 

0<s<t 



(31) 



Proof. We have [Eq. ([23])] H(s) = U (t,O^V(t)U (t,0) 
and [Eq. ©] ll^Wllui = ||J £ W(s)H(s)W(s)Us\\ ui . 
The result follows from the triangle inequality. ■ 

We have presented the bound in the interaction 
picture. Clearly, the same argument applies in the 
Schrodinger picture, where instead one obtains 



l«(*)ll 



1 

< - 

~ t 



ds\\H (s)+V(s)\\ 



(32) 



V. DISTANCE BEFORE AND AFTER PARTIAL 
TRACE 

Since we are interested in the distances between states 
undergoing open system dynamics, we now prove the fol- 
lowing: 

Lemma 2 Let Hs and TIb be finite dimensional Hilbert 
spaces of dimensions d$ and ds, and let A 6 Hs <8> Hb- 
Then for any unitarily invariant norm that is multiplica- 
tive over tensor products the partial trace satisfies the 
following norm inequality 



\\tr B A\\ < 



Ir, 



Mil 



(33) 



where I is the identity operator. 



This result was already known for the trace norm as a 
special case of the contractivity of trace-preserving quan- 
tum operations 

Proof. Consider a unitary irreducible representation 
{UB{g),g € G} of a compact group G on Hb- Then 
it follows from Schur's lemma that the partial trace has 
the following representation [15j |: 

-j-tv B (X) « I B = f[I A ® U B {g)]X[I A ® U B (gYW(g), 

«B JO 

where dfi(g) denotes the left-invariant Haar measure nor- 
malized as J G dfi(g) — 1 and d B = dim(7is). Then 



l|tr B A:|| = 



< 



IKbii jg 

d B 



\tx B (X) ® h 



\\[lA®U B {g)}X[I A ®U B { g y}\\dn{g) 



\I 



B\\ Jg 



\\x\\Mg) 



\x\ 



(35) 



In particular, H/bIIj = d B , \\I B \\ 2 — \/d B , and 
|| Js || oo = 1, and since the trace, Frobenius, and oper- 
ator norms are all multiplicative over tensor products we 
have, specifically: 



\\tr B X\\l < \\X\\ lt 

\\tr B X\\ 2 < y/d^\\X\\ 2l 
||tr B X|l < dsWXlU 



(36) 
(37) 
(38) 



Note that not all unitarily invariant norms are multi- 
plicative over tensor products. For instance, the Ky Fan 
fc-norm ||.||(fc) is the sum of the k largest singular values, 
and is unitarily invariant but it is not multiplicative 
in this way. For example, when d A = ds = d > k > 2 
we have \\I A \\( k ) = ||^s||(fe) = \\Ia ® -^sllrfc) = k. So, 
for X = I A (g> I B we have ||tr B X|| (/c ) = d||J,i||(fc) = dk 



but (di 



L B\\(k), 



X 



(fc) 



(d/k)k — d which gives an 



inequality in the wrong direction. 



VI. DISTANCE BOUND 

We are now ready to prove our main theorem, that pro- 
vides a bound on the distance between states in terms of 
the distance between effective superoperator generators. 
As an immediate corollary we obtain the bound in terms 
of the effective Hamiltonians. 

How much does the deviation between Lit) and L°(t) 
impact the deviation between pit) and p°(t)? We shall 
assume that the desired evolution is unitary, i.e., e tL W = 
U°(t) • r/°(t)t, where U°(t) = e -* tn °W. 

Theorem 1 Consider two evolutions: the "desired" 
unitary evolution p(0) i— > P°{t) = e tL ^ p{0) = 
e - ^*^ W''W0), and the "actual" evolution p(0) i— > 
pit) = e tL Wp(0). Let AL(t) = L(t) - L°(t). Then 

DW), p°it)} < min[l, i( e *H^(*)IU.i _ i)] (39) 



If in addition t\\ AL(t)|| 00j i < 1 then 

DM^A^tllAL^IUi. 



(40) 



Corollary 2 (less tight) For Hamiltonian evolution, 
where L(t) = -§[f2(i),-] and L°(t) = -^[fl°(t) r ], and 
defining Afi(i) = fi(i) - Q°(t): 

D[p(t),p°(t)] < min^ei'll^^-l)] (41) 

< lt\\Anit)\\ x if ^iiao^Hoo < 1. 

n a 

(42) 

Note that in Corollary [2] we use the standard oper- 
ator norm, while in Theorem [1] we use the O-T norm. 
There are two reasons that Corollary[5]is less tight. First, 
it involves an additional use of the triangle inequality, 
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which converts ||[AO(t),p]||i to 2|| AO(i)p||i. Obviously, 
if AQ(t) and p nearly commute then this results in a 
weak bound. Second, even though we do not have an in- 
terpretation of || AL(t)|| 00 i as a function of the singular 
values of AL(t), it is convenient to imagine this to be the 
case; then, for Hamiltonian evolution, the effective su- 
peroperator generator AL(t) has eigenvalues which are 
the eigenvalue (energy) differences of the corresponding 
effective Hamiltonian. Assuming the eigenvalues to be 
positive, their differences are always upper bounded by 
the largest eigenvalue, i.e., || Afi(f)|| 00 . 
Proof of Theorem [TJ Let us define Hcrmitian opera- 
tors 

H° = itL°{t) H = itL{t), (43) 

and consider (superoperator) unitaries generated by 
these operators as a function of a new time parameter 
s (we shall hold t constant): 

dU /ds = -iH Q U , dU/ds = -iHU. (44) 

Then 

U (s) = e- lsU ° U{s) = e-' isW , (45) 
and we can define an interaction picture via 

U{s) = U (s)S(s), (46) 
where the interaction picture "perturbation" is 

v( s ) ee u^(s)(n-n°)u (s) = ml(s)AL(t)u {s). (47) 

Then, using (1) unitary invariance and (2) the definition 
of the || Hoc,! norm [Eq. (fTT|)]: 

D[p(t),p°(t)} = \\\e tL °«Xe- tL °Me tL V -l)p(Q)\\l 
( => ±||(S(1) -2)p(0)||i 
< ilWlJ-JlUx, (48) 

which explains why we introduced S. We can compute 
S using the Dyson series of time-dependent perturbation 
theory: 

dS/ds = -iVS, 

S( S ) = 1 + S rn(s), 
m=l 

S m (s) = / ds-L / ds 2 / ds m V(s 1 )V(s 2 ) • • • V(s m ). 
Jo Jo Jo 

(49) 

Using submultiplicativity and the triangle inequality we 
can then show that (see Appendix IBl for the details): 

||5(a) -JlL^ell^W'lUi-l, (50) 



Thus, finally, we have from Eq. ()48[> : 

D[p(t),p°(t)]<^ tAL ^^-l). (51) 

If additionally t\\ AiH^i < 1 then the inequality e x — 1 < 
(e-l)x yields D[p(t), p°(t)] < tWAl^)]]^. On the other 
hand, note that D[p{t) , p° {t)} = |||(5(1) - T)p(0)\\x < 

1(115(^(0)11! + 11^(0)11!) = 1. - 

Proof of Corollary [2l For Hamiltonian evolution we 
have 

= Bup||[An(t),p]||i 

p 

< sup2||Afi(i)p||i 

p 

< Bup2||An(t)|| 00 ||p||i 

p 

= 2||A0(i)|| oo , (52) 

where in the first inequality we used the triangle inequal- 
ity, and in the second inequality we used Eq. (|6|). ■ 



VII. DISCUSSION 

A. The general bound 

By putting together Eq. (|3"3"|) and Corollaries [TJ [51 we 
can answer the question we posed in the Introduction: 

Theorem 2 Consider a quantum system S coupled to a 
bath B, undergoing evolution under either the "actual" 
joint unitary propagator U(t) = e - "^*^*) or the "de- 
sired" joint unitary propagator U a {t) — e~^ tf2 °'*' ; gen- 
erated respectively by H{t) and H (t). Then the trace 
distance between the actual time evolved system state 
Ps{t) — tr BU(t)p(0)U(ty and the desired one, p^it) = 
tr B U°(t)p(0)U°(ty , satisfies the bound 

D[p s (t),p° s (t)] < min[l,±(et*ll fi W- fi °WIU-l)] (53) 

< min[l, I( e *«ll^ll->-<ll^°ll-» _ 1)]. 

(54) 

where (||-X"||oo) ee \ f Q ds\\X (s) W^. 

This result shows that to minimize the distance be- 
tween the actual and desired evolution it is sufficient to 
minimize the distance in the operator norm between the 
actual and desired effective Hamiltonians, or the differ- 
ence of the time average of the operator norm of the real 
Hamiltonians. Techniques for doing so which explicitly 
use the Magnus expansion or operator norms include dy- 
namical decoupling p^.[l7| and quantum error correction 
for non-Markovian noise[!l, GJ, S El El ■ 
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Illustration using concatenated dynamical 
decoupling 



Acknowledgments 



As an illustration, consider the scenario of a single 
qubit coupled to a bath via a general system-bath Hamil- 
tonian H SB = H s <g> I B + J2 a =x, y , z a <* ® B a (H s is 
the system-only Hamiltonian, a a are the Pauli matri- 
ces, and B a ^ Ib are bath operators), with a bath 
Hamiltonian Hb, and controlled via a system Hamil- 
tonian Hc(t), so that the total Hamiltonian is H(t) = 
He (t) <8> Ib + Hs b + Is ® -^5 ■ Suppose that one wishes 
to preserve the state of the qubit, i.e., we are interested 
in "quantum memory", so that fi°(i) = Is ® It 
was shown in Ref. 17], Eq. (51), that by using concate- 
nated dynamical decoupling (a recursively defined pulse 
sequence [HJ])) and assuming zero- width pulses, the Mag- 
nus expansion yields the following upper bound: 

T\\Q(T) - n°(T)\\ ui /h < JT(/3T/N 1 / 2 ) l ° S4N . (55) 

Here ft°(T) = I s <E> H B , T = Nt (the duration of a 
concatenated pulse sequence with pulse interval t), and 



/3: 



J : 



IB. 



(56) 



are measures of the bath and system-bath coupling 
strength, respectively. Here we shall take the norm in 
these last two definitions as the operator norm. The 
bound ([55]) is valid as long as j3T 1 [l7| • When this is 
the case, the right-hand side of Eq. (|55|) decays to zero su- 
perpolynomially in the number of pulses N . The bound 
(f53|) then yields, for sufficiently large N: 

D[ps(n Ps (T)] < ItfJTWswr** 1} 

< 2JT{[3T/N 1/2 ) los * N , (57) 

which shows that the distance between the actual and 
desired state is maintained arbitrarily well. 



VIII. CONCLUSIONS 

We have presented various bounds on the distance be- 
tween states evolving quantum mechanically, either as 
closed or as open systems. These bounds are summarized 
in Theorem[5J We expect our bounds to be useful in a va- 
riety of quantum computing or control applications. An 
undesirable aspect of Eqs. ([55]) and (|53[) is that the opera- 
tor norm can diverge if the bath spectrum is unbounded, 
as is the case, e.g., for an oscillator bath. A brute force 
solution is the introduction of a high-energy cutoff. How- 
ever, a more elegant solution is to note [24j (lemma 8) that 
every system with energy constraints (such as a bound on 
the average energy), is essentially supported on a finite- 
dimensional Hilbert space. An even more satisfactory 
solution in the unbounded spectrum case is to find a dis- 
tance bound involving correlation functions. This can be 
accomplished by performing a perturbative treatment in 
the system-bath coupling parameter, as is done in the 
standard derivation of quantum master equations [3, HJ ■ 
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APPENDIX A: INTERACTION PICTURE 



We prove that U(t, 0) defined in Eq. (|2"Tj) satisfies 
the Schrodinger equation (|22[) with the interaction pic- 
ture Hamiltonian defined in Eq. (|23[) . This requires 
proof since all the Hamiltonians considered are time- 
dependent, whereas usually one only considers the per- 
turbation to be time-dependent. To see this we differen- 
tiate both sides of Eq. (f21"j) . while making use of Eqs. (|22"|) 
and 



dU(t, 0) 
dt 



U {t,0)U(t,0) 



dt 



dU o {t,0) 



U(t,0)+U o {t, 0) 



dU(t,0) 



dt v ' ' uv ' ' dt 
-iH (t)Uo{t,0)U(t, 0) 
-iU (t,0)H(t)U(t,0) 
-iH {t)U {t,0)U{t, 0) 
-iU (t,0)U (t,0)iV(t)U (t,0)U(t,0) 
-i[H Q (t) + V(t)]U (t,0)U(t,0) 
-iH(t)U(t,0), 



which is the same differential equation as Eq. (fl9| . The 
initial conditions of the equations are also the same 
[U(0, 0) = I], thus Eqs. l[2I ]) -([23" ]) describe the propagator 
generated by H(t). 



APPENDIX B: DYSON EXPANSION BOUND 



We prove Eq. (|5"0| . The proof makes use of (1) the tri- 
angle inequality, (2) submultiplicativity, and (3) unitary 
invariance. 
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\S(s)-Z\L,i 



E s m (s 



rn — 1 

(!) 

< E 

rn— 1 



oo.l 



E 

m— 1 



n 



(1) ~ f s 

< E /„ ds i 

m=l ' 



(I 



(2) ~ 

< E /„ ds i 

m=l ' 



II 



'e/ ^-7 

OO 
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